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Abstract. We prove a law of large numbers for the loss from default and use it for approximating the 
distribution of the loss from default in large, potentially heterogenous portfolios. The density of the limiting 
measure is shown to solve a non-linear SPDE, and the moments of the limiting measure are shown to satisfy 
an infinite system of SDEs. The solution to this system leads to the distribution of the limiting portfolio 
loss, which we propose as an approximation to the loss distribution for a large portfolio. Numerical tests 
illustrate the accuracy of the approximation, and highlight its computational advantages over a direct Monte 
Carlo simulation of the original stochastic system. 



1. Introduction 

Reduced-form point process models of correlated default timing are widely used to measure portfolio 
credit risk and to value securities exposed to correlated default risk. Computing the distribution of the 
loss from default in these models tends to be difficult, however, especially in bottom-up formulations with 
many names. Semi-analytical transform methods have limited scope. Monte Carlo simulation is much more 
broadly applicable but can be slow for large portfolios and longer time horizons. 

This paper develops an approximation to the distribution of the loss from default in large portfolios 
that may have a heterogenous structure. The approximation is valid for a class of reduced-form models in 
which a name defaults at a stochastic intensity that is influenced by an idiosyncratic risk factor process, a 
systematic risk factor process common to all names in the pool, and the portfolio loss rate. It is based on a 
law of large numbers for the portfolio loss rate. The limiting portfolio loss is not deterministic but follows 
a stochastic process driven by the systematic risk factor. This means the exposure to the systematic risk 
cannot be diversified. We show that the density of the limiting measure satisfies a non-linear stochastic partial 
differential equation (SPDE) driven by the systematic risk factor and develop a numerical method to solve 
this equation. The method is based on the observation that certain moments of the limiting measure satisfy 
an infinite system of SDEs. These SDEs are driven by the systematic risk factor; a truncated system can 
be solved using a discretization scheme, for example. The solution to the SDE system leads to the solution 
to the SPDE through an inverse moment problem. It also leads to the distribution of the limiting portfolio 
loss, which we propose as an approximation to the distribution of the loss from default for a large portfolio. 
Estimators of value at risk and other risk measures are immediate from the limiting loss distribution, and 
so are estimators of prices of credit derivatives exposed to the portfolio. 

Numerical tests illustrate the accuracy and computational efficiency of the approximation for large but 
finite portfolios. We find a substantial reduction in computational effort over the alternative of direct Monte 
Carlo simulation of the high-dimensional original stochastic system. The accuracy of the approximation 
mainly depends on the portfolio size N and the sensitivity to the systematic risk factor. For a given 
sensitivity, the accuracy increases with N, as expected. The higher the sensitivity to the systematic risk, 
the higher the variance of the loss distribution and the more accurate is the approximation for fixed N. The 
approximation is remarkably accurate in the tail of the loss distribution, even for relatively small portfolios. 
This renders it particularly suitable for the estimation of portfolio risk measures. 

Large portfolio approximations were first studied by |27j . In Vasicek's static model of a homogenous 
pool, firms default independently of one another conditional on a normally distributed random variable 
representing a systematic risk factor. Because the losses from defaults are conditionally i.i.d., the classical 
law of large numbers ensures the convergence of the portfolio loss rate to its conditional mean, from which 
the limiting loss distribution is immediate. [25] examine alternative distributions of the systematic factor, 
and 123] and [THj study the limiting loss in a heterogenous portfolio. [J] analyze a dynamic extension of 
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Vasicek's homogeneous pool model in which the systematic risk factor follows a Brownian motion. They 
obtain an SPDE driven by that Brownian motion for the density of the limiting measure, which they solve 
using a finite element method. The conditional independence of defaults can also be exploited to analyze the 
tail behavior of the losses in large, not necessarily homogenous, portfolios using large deviations arguments; 
see [H] and dZ]. 

The analysis in this paper differs from that in the aforementioned articles in several important respects. 
We study a class of dynamic point process models of correlated default timing in which a firm defaults at 
a stochastic intensity process. The intensity is influenced by an idiosyncratic risk factor process following a 
square-root diffusion, a systematic risk factor process following a diffusion with arbitrary coefhcient functions, 
and the portfolio loss rate. To address the heterogeneity of a portfolio, the intensity parameters of each name 
are allowed to be different. The choice for dependence of the intensity on the portfolio loss is motivated 
by the empirical results of [2], who find that defaults have a statistically significant feedback effect on the 
surviving firms. The self-exciting behavior of defaults violates the conditional independence property that 
is exploited in the aforementioned articles. It complicates the asymptotic analysis and induces an integral 
term in the drift of the SPDE governing the density of the limiting measure. The exposure to the systematic 
risk leads to the noise term in the SPDE, given by an Ito integral against the Brownian motion driving the 
systematic risk diffusion. The solution to the SPDE governs the distribution of the limiting loss at all future 
horizons, facilitating applications requiring the "loss surface," such as the valuation and hedging of portfolio 
credit derivatives. This dynamic perspective is absent in the static formulations considered in some of the 
aforementioned articles. 

Our law of large numbers significantly extends an earlier result in [15] , which assumes Ornstein-Uhlenbeck 
dynamics for the systematic risk factor and the sensitivity of the intensity to the systematic risk to vanish in 
the limit. The law of large numbers developed here allows for general diffusion dynamics for the systematic 
risk factor. Moreover, the exposure to this diffusion is not required to vanish, generating a much richer, non- 
deterministic limiting behavior governed by an SPDE rather than a PDE. The treatment of these features 
requires additional arguments. The filtered martingale problem is used to identify the limit and prove the 
law of large numbers. A major difficulty in the identification of the limit is the solution of a coupled system of 
SDEs, which we address using fixed-point arguments. In contrast to [15], this system does not easily decouple 
in the more general setting considered here. The analysis of the fixed-point arguments is complicated due 
to the square-root singularity. 

[B] prove a LLN for a related system, taking an intensity as a deterministic function of an idiosyncratic risk 
factor, a systematic risk factor, and the percentage portfolio loss. The risk factors follow diffusion processes 
whose coefficients may depend on the portfolio loss. In that formulation, the impact of a default on the 
dynamics of the surviving firms is permanent. In our work, an intensity depends on the path of the portfolio 
loss. Therefore, the impact of a default on the surviving firms may be transient and fade away with time. 
There is a recovery effect. Other interacting particle systems with permanent default impact are analyzed 
by [8] and [7], who take an intensity as a deterministic function of the portfolio loss rate. In a model with 
local interaction, |16) take the intensity of a name as a deterministic function of the state of the names in a 
specified neighborhood of that name. These papers prove LLNs for the portfolio loss and develop Gaussian 
approximations to the portfolio loss distribution based on central limit theorems. The interacting particle 
system which we propose and study includes firm-specific sources of default risk and addresses an additional 
source of default clustering, namely the exposure of a firm to a systematic risk factor process. This exposure 
generates a random limiting behavior. 

There are several other related articles. [10] develop large portfolio approximations based on a law of large 
numbers and a central limit theorem in a stochastic network setting, in which firms default independently 
of one another conditional on the realization of a systematic factor governed by a finite state Markov chain. 
[26] examine large portfolio asymptotics for utility indifference valuation of securities exposed to the losses 
in the pool. Our formulation addresses the dependence of the intensity on a systematic diffusion factor and 
the portfolio loss. It allows for self-exciting effects that violate the conditional independence assumption. 
Gaussian and large deviation approximations to the distribution of portfolio losses for an affine point process 
system with features similar to those of our, not necessarily affine, system arc provided by [29j . Their 
asymptotic analysis is based on a "large horizon" rather than a "large portfolio" regime considered here and 
in the aforementioned articles. The scope of their approximations differs from that of ours. 
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The rest of the paper is organized as follows. Section [5] introduces our model of correlated default timing 
in a portfolio of names. Section [3] states our main result, Theorem 13. 1[ a law of large numbers for the loss 
rate in the pool. Section 3] provides further insights into the limiting behavior of the loss for the special case 
of a homogenous pool. Section [5] develops and implements a method for the numerical solution of the SPDE. 
Numerical results illustrate the method and demonstrate the accuracy of the approximation of the portfolio 
loss by the limiting loss. Section [5] discusses the extension of our results to more general intensity dynamics. 
Sections [7] and E] are devoted to the proof of Theorem 13. II Appendices provide auxiliary results. 

2. Model and Assumptions 

We provide a dynamic point process model of correlated default timing in a portfolio of names. We assume 
that (fi,^,P) is an underlying probability space on which all random variables are defined. Let {VF"}ragN 
be a countable collection of standard Brownian motions. Let {e„}„gN be an i.i.d. collection of standard 
exponential random variables. Finally, let V he a, standard Brownian motion which is independent of the 
W^'s and e„'s. Each W" will represent a source of risk which is idiosyncratic to a specific name. Each e„ 
will represent a normalized default time for a specific name. The process V will drive a systematic risk factor 
process to which all names are exposed. Define V = and = J^^-^ , the completed subfiltrations of 
^ generated by V and {V, W^), respectively. We have V = IJtGK 

Fix N G N, n E {1,2, . . . , N} and consider the following system: 

-ajv,„(Af ^" - ~XN,n)dt + aN^n\/)^dWr + l3'^,„dLf + /3^_„Af "dX* t > 

Ao,Af,Tl 

ba{Xt)dt + <Jo{Xt)dVt t>0 

Here, x is the indicator function. The initial condition Xo of X is fixed. The description of is equivalent 
to a more standard construction. In particular, define 

(2) r^^" ='inf |t>0:^* Af '"ds > e, 

Then X[c„,oo)(X*=o Af '"ds) = X{r" "<i} and consequently 

1 ^ 

(3) ]^H^{^"'"<*}- 

The process represents the loss rate in a portfolio of N names, assuming a loss given default of one 
unit. The process A^'" represents the intensity, or conditional event rate, of the n-th name in the pool. 
More precisely, A^'" is the density of the Doob-Meyer compensator to the default indicator X{T"."<t}; see 
(jni). Lemma 3.1 in [TS] implies that the A^'"'s are indeed nonnegative. 

The jump-diffusion intensity model ([Ij is empirically motivated. It addresses several channels of default 
clustering. An intensity is driven by an idiosyncratic source of risk represented by a Brownian motion 
VF", and a source of systematic risk common to all firms-the diffusion process X. Movements in X cause 
correlated changes in firms' intensities and thus provide a channel for default clustering emphasized by [9] 
for corporate defaults in the U.S. The sensitivity of A^'" to changes in X is measured by the parameter 
Pf; n ^ The second channel for default clustering is modeled through the feedback ("contagion") term 
ndL^ ■ A default causes a jump of size „ in the intensity A^'", where „ € IR+ = [0, oo). Due to 
the mean- reversion of X^-^ , the impact of a default fades away with time, exponentially with rate o. m ^ . 
[2] have found self-exciting effects of this type to be an important channel for the clustering of defaults in 
the U.S., over and above any clustering caused by the exposure of firms to systematic risk factors. 
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(1) 



dXt = 
Xo = 



L 



N 



We allow for a heterogeneous pool; the intensity dynamics of each name can be different. We capture 
these different dynamics by defining the "types" 

,711 ^N.n: I^N.n^ /^A^,n)' 

the p^'"'s take values in parameter space V =^ x M. In order to expect regular macroscopic behavior 
of as N oo, the p^'"'s and the Ao,7v,n's should have enough regularity as — ?• oo. For each N G N, 
define 

N N 
71—1 n — 1 

these are elements of !3^(V) and ,iF(M+) respectiveljQ- 

We require three main conditions. These conditions are in force throughout the paper, even though this 
may not always be stated explicitly. Firstly, we assume that the types of ^ and the initial distributions 
(the Ao, Aran's) are sufficiently regular. 

Condition 2.1. tt limjv^oo and Ko ="'^ lim^v^oo exist [in ,^{V) and ^(R+), respectively). 

We also require that the tt^'s and A^'s all (uniformly in N) have compact support. Wc could relax this 
requirement, at the cost of a much more careful error analysis. 

Condition 2.2. There is a K > such that the aN,n 's, XN,n 's, crN,n 's, „| 's, „| 's, and Xo.w.n 's are 
all bounded by K for all N £N and n€ {1,2, ... ,N}. 

Regarding the systematic risk process X, we assume 

Condition 2.3. The functions bo and cfq that determine the systematic risk process X are such that the 
corresponding SDE has a unique strong solution. Moreover, there is a T > and a function u{x) such that 
aoix)u{x) = —bo{x) and satisfying 

(5) E[e^/o"l"(^=)lH 



< oo. 



The Novikov condition ^ may not be necessary. Lemma 18.21 is the key step for the proof of a law of 
large numbers for the loss rate in the system ([1]), which is stated as Theorem 13. II Its proof is based on 
a fixed point argument and uses Girsanov's theorem; this is where ([5]) is required. 

Our basic formulation significantly extends that of [15]. First, we allow the systematic risk X to follow a 
general diffusion process with coefficients satisfying Condition 12.31 rather than a simple Ornstein-Uhlenbeck 
process. Second, we no longer require the exposure to the systematic risk, ^, to vanish in the limit as 
N — >■ oo. This implies a richer, non-deterministic limiting behavior. The analysis of this behavior is more 
challenging and requires new arguments. 

Section [S] discusses further extensions of our basic formulation, including stochastic position losses and 
more general intensity dynamics. 

3. Law of Large Numbers 

We develop a law of large numbers for the portfolio loss rate in the system ([T]). To this end, we need 
to understand a system which contains a bit more information than the loss rate . For each N E N and 
n€ {1,2, ...,iV}, define 

(6) Mf X[o,c„) Af ^"d.) = X{.«,.>n 

(where r^'" is as in In other words, ivif^'" = 1 if and only if the n-th name is still alive at time t; 

otherwise Mj ~ 0. Thus m^^" is nonincreasing and right-continuous. It is easy to see that 



^As usual, if B is a topological space, ,0^(E) is the collection of Borel probability measures on E. 
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is a martingale. Define V'^^^V xR+. For each N eN, define pf (p^^", Af '") for all n e {1,2, . . . , N} 
and t > 1. For each t > 0, define 

1 ^ 

n=l 

in other words we keep track of the empirical distribution of the type and intensity for those assets which 
are still "alive" . We note that 

Lf = \~^if{v), t>o. 

We want to understand the dynamics of /i^ for large N (this will then imply the "typical" behavior 
for L^). To understand what our main result is, let's first set up a topological framework to understand 
convergence of /j,^ . Let E be the collection of sub-probability measures (i.e., defective probability measures) 
on V; i.e., E consists of those Borel measures v onV such that J^('P) < 1. We can topologize E in the usual 
way (by projecting onto the one-point compactification of V; see [24l Ch. 9.5]). In particular, fix a point * 
that is not in V and define =^ 'P U {★}. Give the standard topology; open sets are those which are 
open subsets of V (with its original topology) or complements in of closed subsets of V (again, in the 
original topology of V). Define a bijection t from E to ^('P^) (the collection of Borel probability measures 
on ■p+) by setting 

(ii.)(A) v{A f^V) + {l- v{V)) 5M) 

for all A e 3§{V^). We can define the Skorohod topology on ,'3^{V^), and define a corresponding metric on 
E by requiring t to be an isometry. This makes E a Polish space. Thus, //^ is an elemen10 of -D_e[0, oo). 

The main result of this paper is Theorem 13.11 essentially a law of large numbers for &fi N ^ oo. For 
p = (p. A) where p — (a. A, cr, (3'-' , (3^) G V and / e C°°{V), define the operators 

iCm) = ^y^Q^iP) - «(A - A)^(p) - A/(p) 

(/:2/)(p) = |{(p) 

(7) 

{Cym = /3^A6o(x)-((p) + {/3S^^X'al{x)-^{p) 



(£^/)(p)^/3^Aao(x)|{(p). 



Also define 



The generator Ci corresponds to the diffusive part of the intensity with killing rate A, and £2 is the macro- 
scopic effect of contagion on the surviving intensities at any given time. Operators £3 and £| are related to 
the exogenous systematic risk X. 

Theorem 3.1. We have that fi^ converges in distribution to j2. in De[0,T]. The evolution of p,. is given 
by the measure evolution equation 



+ (^cf'f,fit)^dVu yf ec°^\r) a.s. 



dt 



E 

Suppose there is a solution of the nonlinear SPDE 

dv{t, p) = S^Clv{t, p) + Cl'^^'vit, p) + (^^^ ^ Q{p')v{t, p')dp'^ C;v{t, p) \ dt 
+ Cl'^'v{t,p)dVt, t>0, peV 



■^If 5 is a Polish space, then _Ds[0, 00) is the coUeetion of maps from [0, 00) into S which are right-continuous and which 
have left-hand limits. The space D^^O, 00) can be topologized by the Skorohod metric, which we will denote by dg; see I12| . 
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where C* denote adjoint operators, with initial condition 

lim v{t, p)dp = TT X Aq. 

Then 

P-t = v{t, p)(ip. 

Remark 3.2. The SPDE should he supplied with appropriate boundary conditions. In Section \4.2\ below, 
we will justify the conditions 

v{t, A = 0, p) = v{t, A = oo, p) = 0. 

The proof of Thcorcm l3.1l is given in Sections [7] and [51 Lemma |8 . 41 provides an alternative characterization 
of the hmit /l. AuxiUary resuhs are given in the appendices. 

Remark 3.3. Equation is a stochastic partial integral differential equation (SPIDE) in the half line that 
degenerates at the boundary A = 0. Due to Lemma \8.Si below, (0) can be viewed as a linear SPDE in the half 
line that degenerates at the boundary. Indeed, by Lemma \8.2l there is a unique pair {{Q(t), Xt{p)) '■ t > 0} 
taking values in x satisfying the coupled system ll28\) -[2 ^) . Noting that Q{t) = {Q,flt)E, we see that 
the SPDE (0) can be written as 

(9) dv{t,p) = [civit,p) + C;''''v{t,p) + Q{t)C*Mt,p)^dt + Cl''''vit,p)dVt, t>0, p€V. 

Notice also that by Remark \4-.1\ Q{t) is bounded for all t £ IR+. Linear SPDEs in the half line that degenerate 
at the boundary are treated, under alternative assumptions, by [21) . [20j . |22| and |19| . 

4. Homogeneous Pool 

We develop further insights into the SPDE governing the limit density in the case that the portfolio is 
homogenous. Let p (p, A) where p = {a,\,a, (3^ , (3^) G V. For a homogenous pool, p^'" = p for all 
iV € N and n £ {1,2,..., N}. In this case, we write v{t, A) for the solution of the SPDE dS]), suppressing 
the dependence on the fixed p. 

4.1. Limiting Portfolio Loss. The SPDE takes the form 

dv{t, A) = \civ{t, A) + Cl'^'v(t, \)+(3^ ( [ \v{t, \)d\ \ Clv{t, A) j dt + C*^^'v{t, X)dVt, t,X>0 



(10) u(0,A) Ao(A), 



v{t, 0) = lim v{t, A) = 
Ayoo 

where the adjoint operators are given by 

32 



C*2vit,\) 

(11) 



^*Mt^ ^) = ^ A)^ + ^ («(A - ~X)vit, A)) - Xvit, A) 

dv{t,X) 



dX 

32 



= ^ (^lif3'rX'aUxHt,X)^ ^ (/3«A6o(x)^;(^, A)) 



Clv{t,X) = -(3'''ao{x)^{Xv{t,X)) . 
Define the limiting portfolio loss L by 

(12) Lt 1 - / v{t, X)dX, t > 0; 



this is a random quantity since v{t, A) depends on the systematic risk Xf . For large N, Theorem [3TT] suggests 
the "large-portfolio approximation" 

(13) if w it, t> 0. 
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4.2. Justification of Boundary Conditions. Since v{t, X)dX < 1 and v{t, A) > for alH > 0, A e M+, 

limA-s-oo vit, A) = 0. 



The boundary condition of v{t, A = 0) is implied by the intensities Af^'" being positive almost surely. For 
the deterministic case of ~ 0, it is sufficient to stipulate Feller's condition of aX > ^cr^, see [13]. For the 
case of (3^ > 0, let us assume that v{t,X = 0) < oo. Then, Feller's condition is again sufficient to imply 
v{t, A = 0) = 0. Let the flux be f{t) where f(t) =^ ^ v(t, X)dX = - Xv{t, X)dX. This follows from 
the fact that the Af stay non-negative almost surely by Lemma 3.1 in [15] and therefore, according to 
the empirical measure fi^ , only leave [0,oo) by defaulting. In the asymptotic case described by the SPDE, 
defaults occur via the sink term —Aw. Then, probability mass only leaves [0, oo) via the sink term and does 
not flow across the boundaries at A = and A = oo. Integrating pH)) over and using the aforementioned 
flux condition along with the boundary condition at A = oo, we have that 

= aXv(t,X = 0)- -a^v{t,X^0)+l3^{ / Xv{t, X)dX)v{t, X ^ 0), 
2 Jo 

which is only satisfied by v{t, A = 0) = 0. 

This discussion provides a justification for the choice of the boundary conditions for a homogeneous pool. 

The treatment of a heterogeneous pool is analogous. 

4.3. Alternative Representation of Limiting Loss. As we shall see below, if the density v{t, A) has 
sufficiently fast decay at A = oo and at A = 0, then one can justify an alternative representation of the 
limiting loss p2)) . By integration by parts, we have 



1 - 

1 - 



v{t,X)dX 



v{0,X)dX+ I [a{X~X)v{s,X)~^'^Xbo{Xs)v{s,X)]lJ^ds 
Jo 

^<J^Xvx{s, A) + ^<J^v{s, A) + i X^alvxis, A) + {p'^y X4v{s, A) 



ds 



A=0 



(14) 



Xv{s,X)dXds - 13^ 



'0 JO 

Now, observe that 



This implies that 
(15) 



(t,u(s,-)) [w(s,A)]^^;^ ds 



[/?^Aao(A,)i;(s,A)]J:;;°dl/. 



N N 
n=l 



N ^ Po 

n=l 



>>0} 



= 1. 



u(0,A)dA = 1. 



Recall the boundary conditions v{t, 0) = limxy^ao v{t, A) = 0. Then, if we assume that for any t € M+, v{t, •) 
and v\{t, •) decay fast at infinity, in the sense that lim^^oo A^UA(i, A) = lim^^oo Xv{t, A) = 0, the integration 
by parts formula (fT4]) and equation (fTSj) imply 



(16) 



1 



v{t,X)dX 



Xv{s, X)dXds 



A particularly interesting consequence of (|16p is summarized in the following remark. 

Remark 4.1. Relation il6]) implies that the integral term in ilO\) is bounded. Indeed, the term is Q(t) 
l3^Xv{t,X)dX. Then, (01) can he rewritten as 



(17) 



1 

u(i,A)dA + -^ 



Q{s)ds = 1 



Since v{t, X)dX < 1, relation Ji?] ) implies that for any t G M_|_ we have Q{t) < oo. It is easy to see that 
the corresponding conclusion also holds for the heterogeneous pool due to Condition \2.2[ 
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Moreover, (fT5|) implies for the rate of change 



oo 



Lt^~ vt{t,X)dX^ / Xv{t,X)dX^ / A/2t(dA). 
JO ^0 Jo 



Finally, we mention that one can view {Lt,v{t, A)) as a pair satisfying 

dv{t,X) = ^^Clv{t,X)+Cl'^*v{t,X)^dt + p'^'C2v{t,X)dLt+Cl'^'v{t,X)dVt t,X>0 



u(0,A) = Ao(A) 

Lt = 1-/ v{t,X)dX^ / / Xv{s,X)dXds. 
Jo Jo Jo 

5. Numerical method and results 

Wc develop and implement a method for the numerical solution of the SPDE pop governing the limiting 
portfolio loss p2|) in a homogenous pool. We obtain the distribution of the limiting loss Lt, which we propose 
as an approximation to the loss Lf^ when N is large. Numerical results illustrate the method, as well as the 
accuracy and efficiency of the approximation. 

5.1. Numerical Approaches in the Deterministic Case. In the case of (3^ = 0, pO)) becomes a deter- 
ministic, quasi-linear PDE. For completeness, we outline some numerical methods for this case. 

5.1.1. No Feedback. If in addition (5^ = 0, the default times r^'" are independent and we have the analytic 
solution Lt = l- exp{A{t) + B(t)A)Ao(A)dA where A(t) and B{t) satisfy 

B{t) = ^ (^a + jta.nh(^-^jt + ci 

Ait) = -C2d2t H log f— 

7 \ cosh(cij 

where 7 = \/kM~2ct2, ci = tanh~"'(^), C2 = di ~ and d2 = — aA. 

5.1.2. Finite Difference. For the deterministic case where (i'-^ > 0, a finite difference scheme can efficiently 
solve the PDE. We devise a scheme which is implicit in the differential operators and explicit in the integral 
operator. A predictor-corrector iteration is employed to increase accuracy for the integral term. This finite 
difference scheme is second-order accurate. Formally, 



A 



^-^^^^ = C[\{v,+,+v,)]+I[\{v^^i^+v,)]Cl[\{v,+,+v,)l 
where C = C\+Cl and T[v{t, A)] Xv{t, X)dX. 



5.2. Method of Moments. Wc provide a method for the numerical solution of the SPDE (|10p that applies 
in the case that (3^ > 0. Suppose that the boundary conditions for the SPDE are v{t, A = 0) = and 
lim>,_j.oo 'v{t, A) = 0, as justified in Section 14.21 above. (Note that the latter boundary condition also implies 
limA-).oo v\{t, A) = 0.) Furthermore, suppose that for each fc S N, X''v{t, A) is integrable on R+, almost surely. 
Then, the moments Uk{t) ~ X'^v{t, X)dX exist almost surely. A sufficient condition is that the solution 
v{t,X) decays exponentially in A; that is, there exist constants Ci,C2 > such that v{t,X) < Cie~'~^^^ 
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almost surely for t > 0, \ E R+. (We note that it was shown in Remark 14.11 that uq and ui exist.) Then, 
the moments follow the SDE system 

dukit) = {ukit){ -ak + l3^boiXt)k + 0.5ip^ fal{Xt)kik ~ 1)) 

(^g) + itfc_i(i)(0.5a2fc(fc - 1) + aXk + (3^ kui{t)) - Uk+i{t)]dt + p'^ cTo{Xt)kuk{t)dVt, 

uk{t = 0) = / A'^Ao(A)(iA. 



^0 

To find ui;{t), multiply (jlOp by A*^ and integrate by parts over [0,oo). Also, use the boundary conditions at 
A = and A = oo. Note that the limiting loss Lt ~ 1 — UQ{t). 

Moment methods have previously been applied to deterministic PDEs such as the Boltzmann equation, 
see [28] ■ for example. The moment SDE system in our case is not closed since the fc-th equation introduces 
the {k + l)-th moment. So, in practice one must perform a truncation at some level k = K where we let 
uk+i = UK (that is, we use the first A' + 1 moments). The incremental error in the deterministic case 
= 0) resulting from not including the fc-th moment is of the order The error in the stochastic case 
{(3^ > 0) is harder to analyze. Numerical results, reported below, will show that the convergence in terms 
of the number of moments is very rapid. 

We remark that /3'~^ plays a pivotal role in the truncated version of system ([T5|) . If ~ and appropriate 
choices are made for the coefficient functions 6o(') ^nd (7o(-) of the systematic risk, the truncated system 
satisfies the global Lipschitz condition and we have the standard existence and uniqueness results. If 13^ > 0, 
the truncated system is only locally Lipschitz and therefore there exists a unique solution up to a stopping 
time C(uo(t = 0),...,UK{t = 0),w) : R^+^ xfl — > [0,oo). 

In addition, the moments can be inverted to yield v{t,X). However, although a distribution uniquely 
determines its moments, the converse may not be true. See [1] for some numerical methods for moment 
inversion. 

There is an alternative approach to viewing the moment system. The system (|18p is driven by a single 
diffusion, which suggests there should be a canonical form where only one SDE has a diffusion term and the 
other SDEs only have drift terms. This is advantageous for numerical solution. In the original system (|18p . 
there is potential for bias when integrating the UkdVt term. In the canonical form of the system, these terms 
are removed and the only source of this type of bias is in the simulation of the process X . The process X 
can, however, often be simulated exactly; see [3, [5], and [13] for exact simulation methods for SDEs. To 
find the canonical form of the moment system, define the new variables rik{t) = Xt — \og{uk{t)) for fc > 1 
and > 0. We get 

dXt = bo{Xt)dt + aoiXt)dVt, 

duoit) = -e'5'(^''"i(*))di, 

dVkit) = {bo{Xt) + hp''al{Xt)^^[^a + p''boiXt) + ^{p'^fal{Xt){k-l) 



ig(fc+l)/3''(Xt-,,fc + i)-fc/3S(Xt-,,fc)-| 



}dt. 



5.3. Behavior of Limiting Loss Distribution. We provide some numerical results. Here and below, we 
choose the common risk process to be the CIR process dXt = k{9 — Xt)dt + eyOCtdVt where k = 4, = .5, 
e = .5, and Xq = .5. We choose the initial condition Ao(A) = S{X — Xq). The moment system ^TE\\ is solved 
using an Euler scheme with time-step of .01. Figure [1] shows the rapid convergence of the moment system 
solution. Even using as few as six moments {K = 5), one can achieve a very accurate distribution for the 
limiting loss Lf. Also, it is noteworthy that in the deterministic case of fS^ = 0, the moment method is faster 
for the same accuracy than the finite difference approach outlined in Section fS. 1.2 1 

We report some salient features of the limiting loss distribution. Figure shows the effect of the feedback 
sensitivity parameter on the distribution of Lt- As /3'~' increases, the mean of the losses increases and a 
heavy tail develops on the right (indicating a greater probability of extreme losses). Larger (3^ also causes 
a wider or more spread-out distribution, indicating a higher variance. An important ramification is that 
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Figure 1 . Comparison of distribution of limiting portfolio loss Lt for different truncation 
levels K at t = 1. The parameter case is cr = .9, a = 4, A = .2, Ao = .2, = 2, and 
= 2. 50, 000 Monte Carlo trials were used. 

greater connectivity between firms, modeled here through a nonlinear term, can increase the volatility of 
their ensemble behavior. Similarly, increasing the systematic risk sensitivity parameter causes heavy tails 
on the right, see Figure [H A hypothesis is that the initial losses can be sparked by the systematic risk factor 
(whose influence is determined by the parameter (3^) and then are later magnified by the contagion risk 
factor (determined by the parameter Z?*^). This has significant economic implications for the spread of risk 
through macro credit markets. The joint presence and interaction of systematic and contagion risk greatly 
magnifies the likelihood of extreme default events. Another interesting observation is that the solution can 
demonstrate non-monotonic behavior in time; for example, the distribution of losses may widen and then 
later tighten. Figure S] shows the evolution of the distribution of Lt over time i, demonstrating this non- 
monotonicity. It is noteworthy that our numerical method yields the limiting loss distribution for all horizons 
simultaneously; this is useful for some applications, including the valuation of portfolio credit derivatives 
(such as index and tranche swaps) and the analysis of portfolio risk measures such as value at risk. 

5.4. Accuracy of Large-Portfolio Approximation. We analyze the accuracy of the approximation (|13p . 
To this end, we estimate the distribution of the loss in a pool of N names by Monte Carlo simulation 
of the default times r^'". The simulation uses a time-scaling method, which is based on a discretization 
of the intensities A^'" on a time grid {tjj^^i where tj ~ jA for some A > 0. Intensities arc simulated on 
using a truncated Euler scheme so that they remain nonnegative. Firm defaults (as well as jumps 
in intensity due to defaults) occur at the grid points according to a discretized version of ([2]): 

• Generate independent e„ ^ Exp(l) for each firm n = 1,2, . . . , N, 



• end. 

The increments Xtj_^-^ — Xt- can be simulated using an Euler scheme similar to the one described above or 
by an exact scheme. For the results presented here, we again choose A" to be a CIR process with the same 
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.^■") V2 V^(0, 1) + /3^A^^"(A,^.,, - A,J] , 



) > e„, then t^'" = tj+i. 



60 




0.55 



Limiting Portfolio Loss 

Figure 2. Comparison of distribution of limiting portfolio loss Lt for different values of 
the contagion sensitivity P'^ at i = 1. The parameter case is cr = .9, a = 4, A = .2, Ao = .2, 
and (5^ = 2. 15, 000 Monte Carlo trials and 16 moments {K = 15) were used. 

parameters as stated earlier. The CIR process X is simulated using an Euler scheme (truncated at zero as 
shown above). Wc choose a time-step of A = .01. 

Convergence of the distribution of to that of the limiting loss Lt tends be more rapid when (3^ is larger. 
When the variance of the losses is very small (i.e., the limiting losses arc close to a deterministic solution), 
the convergence rate is slower. Convergence is most rapid in the tails of the loss distribution. Figures [5] and |6] 
show the convergence of the distribution of to that of Lt for two different parameter cases. Convergence 
is relatively slow, but does indicate that the asymptotic solution is applicable for a portfolio consisting of 
several thousand names, a portfolio size not unusual in practice. Figure [7] compares the value at risk (VaR) 
of L^ and Lt at the 95 and 99 percent levels. The limiting VaR is surprisingly accurate even for smaller 
portfolios. This indicates the usefulness of the approximation for risk management purposes, which focuses 
on the tail of the loss distribution. 

We evaluate the computational efficiency of the approximation. Computations arc performed for the case 
of (T = .9, a = 4, A = .2, P'^ = 2,(3^ = 3, Aq = .2, and t = 1. A discrete time-step of .01 is used. 1000 Monte 
Carlo samples are produced. In general, the computational effort will be of an order greater for the 
generation of L^ , where N is the number of firms in the pool and K is the truncation level of the moment 
method. The computations were done in the Matlab environment on a computer running Mac OS X with a 
2.7 GHz dual-core processor. 



Method of Moments with K = 200 


6.25 seconds 


Finite System with N = 500 


13.58 seconds 


Finite System with N = 1000 


22.12 seconds 


Finite System with N 5000 


236.24 seconds 


Finite System with N = 10000 


441.80 seconds 


Finite System with N = 25000 


1048.80 seconds 



6. Extensions 
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Figure 3. Comparison of distribution of limiting portfolio loss Lt for different values of 
the systematic risk sensitivity ai t = 1. The parameter case is tr = .9, a = 4, A = .2, 
Ao = .2, and = 2. 15000 Monte Carlo trials and 16 moments {K = 15) were used. 



6.1. Extending the model. We can extend the system ([Ij to the case of more general coefBcient functions 
for the intensity as well as stochastic position losses. Fix € N and n £ {1, 2, . . . , N}. Let be 
a family of positive and bounded random variables. The variable d^'" represents the loss at default of the 
n-th name in a pool of size N. Consider the following system: 

dAf = biX^'")dt + a{X^'")dWr + l3%J{Xt)dLf + /3^_„7(Af t > 

\N,n _ 

dXt ^ ba{Xf)dt + ao{Xt)dVt t > 

Xq — Xo 

N 



(19) 



tN 1 \ " gN,n 



N 

n=l 



where the coefficient functions b{-), cr(-), 9{-) and 7(-) satisfy suitable regularity conditions guaranteeing the 
stence of a unique no 
Define the operators 



existence of a unique nonnegative solution 



(/:2/)(p) = |{(p) 

(cum - /3'i{\)boix)^ip) + ^{i3'ri'{x)<jUx)^{p) 

(£^/)(p) = /3S(A)ao(x)|{(p) 
Q(p) = 
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Figure 4. Evolution of distribution of limiting portfolio loss Lt over time t. The parameter 
case is a = .9, a = 4, A = .2, Ao = .2, = 4, and 13^ = 2. 15, 000 Monte Carlo trials and 
16 moments {K — 15) were used. 

Then, following the arguments applied to the system ([1]), we can show that the SPDE governing the limiting 
density takes the following form: 

dv{t, p) = f^civit, p) + p) + e{Xt) ^^^^ Qip'Mt, p')rfp') c;v{t, p)| dt + ci^'^'vit, p)dVt, 

for t > and peV, where p = (/3^, /J-^, e P R+ x M x R and P '^=^V x M+. 

Naturally, one expects that, for large N and for every t > 0, i ~ v{t, p, X)dX where £ = 

6.2. Extending the Moment Method. The moment SDK system can be extended to the case of gen- 
eral coefficient functions as well as the non-homogeneous parameter case. By Theorem 13.11 in the non- 
homogeneous case the SPDE takes the form 

dv{t, p) = \^^a^{Xv{t, p))aa + a((A - ~\)v{t, p))x - \v{t, p) - /3^6o(^t)(Aw(i, p))a 
+ ^(/3^)'^o(^t)(A'?^(t,p))AA-/?^ / Xv{t,p)dXdp{v{t,p))x}dt 

- f o^{Xt){Xv{t,p))xdVt, 
A = 0, p) = A = cx), p) = 0, 

u(t = 0,p) - h{p). 
The SDE moment system for the non-homogeneous case follows: 

p) = p)\-ak + fh^{Xt)k + \^{p^fal{Xt)k{k - l)]dt - p)dt 

+ Ufc_i(t, p)[0.5CT^fc(fc- 1) +aAfc + /3'^fc / UY{t,p)dp\dt + P^aQ{Xt)kuk{t,p)dVt, 

/>oo 

ufc(t = 0,p) = / X''h{p)dX, 
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Portfolio Loss 

Figure 5. Comparison of distributions of limiting portfolio loss Lt (using 16 moments) 
and portfolio loss of finite system for different N a.t t = 1. 25,000 Monte Carlo trials 
were used for the finite system and 100, 000 for the asymptotic solution. The parameter 
case is CT = .9, a = 4, A = .2, Aq = .2, /3'^ = 2, and /J'^ = 3. 

where Uk{t, p) = X''uk{t, p)dX. The SDE moment system is coupled across the parameter space V. For 
numerical implementation, the parameter space V must be discretized. 

For the case of intensity processes with general coefficients b{-), ct(-), and 7(-), we have an SPDE of the 
form stated in Section l6Jl A general moment method can be applied to this class of SPDEs, provided that 
we prescribe b{-), cr(-), and 7(-) such that the processes A^'" stay positive almost surely. Also, assume 6(-), 
(t(-), and 7(-) are analytic on R+. Then, the generalized moment 

uuMMM.MMMit) = I ^'^(A)'"al^-(^)''^5X^7(A)'^-^^«(i,A)dA 
solves a moment system similar (albeit more complicated) to (|18p. 

7. Tightness and Identification of the Limit 
We start by discussing relative compactness of the sequence {/i^jArgN. 
Lemma 7.1. The sequence {fi^jNefi is relatively compact in D^p, oo). 

Proof. The proof of this lemma follows exactly as in [TS]. Thus, the details are omitted. □ 

Next, we want to use the martingale problem (see [H]) to identify the limit of /i^'s. For every / S C°°{V) 
and II IE E, define 

=' / fip)Kdp)- 

Let S be the collection of elements $ in i3(R x ^{V)) of the form 

(20) ^"(2;,^) = fi{x)f2 {{fi,n)E, (./2,A*)£; ■ • • {fM^t^E) 

for some A I e N. some ipi e C°°(R), some e C°°(R^O and some {/mji^i. Then S separates ^(K x V) 
[12] . It thus suffices to show convergence of the martingale problem for functions of the form ((20|) . 
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Figure 6. Comparison of distributions of limiting portfolio loss Lt (using 16 moments) 
and portfolio loss of finite system for different N a.t t = 1. 25,000 Monte Carlo trials 
were used for the finite system and 100, 000 for the asymptotic solution. The parameter 
case is CT = .9, a = 4, A = .2, Aq = .2, /3'^ = 4, and 13^ = 8. 

Let's fix / e C°°{V) and understand exactly what happens to (/, M^)^; when one of the firms defaults. 
Suppose that the n-th firm defaults at time t and that none of the other names defaults at time t (defaults 
occur simultaneously with probability zero). Then 

l<n'<N \ / 
n 

l<n'<N 
n' 

Note furthermore that the default at time t means that J^*^^ Xf'"ds = c„, so M^^'" = 0. Hence 

(21) (/,A*f)B-(/,Mf-)^ = J-^.JO 

where 

U\ '^^^ ^ i f f ^N,n' \N,n' , I^N,n\ r f „N ,71' \ N .n' \\ ^ N ,n' 1 r ( N .11 ,n\ 

= ^ j/l P ' +^j-/(^p',A, jjM, -_/(^p ,A, j 

for ah t > 0, e N and n e {1, 2, . . . , N}. 
For convenience, let's define 

for all / G C^(M). This is the generator of the systematic risk. 

We now identify the limiting martingale problem for fi^ . For p = (p. A) where p — {a,X,a,P^,P^) G V 
and / e C°°{V), recall the definitions of the operators in ([7]). 

For $ G 5 of the form ([20|) we define the following operators. 
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Figure 7. Comparison of 95 and 99 percent VaR of distributions of limiting portfolio loss 
Lt (using 16 moments) and portfolio loss of finite system for different N and several 
horizons t. The parameter case is cr = .9, a = 4, A = .2, Ao = .2, (3^ = 2, and /S*^ = 4. 
50,000 Monte Carlo samples are used for the finite system and 250,000 for the asymptotic 
solution. 



M 



X] i{fl,f^)E^ih^f^)E---ifM,fJ-)E) {{Mx)Clfm, Ai)i3} 



and 



(i3$)(a;,/^) =^ ao{x)^^{x)^p2{{fl,^J)E,{^2,^J)E■■■{hI■,^AE) 

(22) +¥^1(2:)^ ^((/i,a*)^,(/2,a*)£...(/m,m)^){(/:^/,„,m)^} 

m=l 

Moreover, we define the following processes 
(23) 

TV 

(i^^"$)(X„ /.f ) '^=^f ^i(X,) ^ Af {^2 ((/l, A^f )^ + Jknit)^ (/2, A^f )^ + (t) . . . (/m, /if )^ + Jl,%t) 



n=l 



m—1 



16 



■(pf^")Mf''"dw^; 



and 
(24) 

m— 1 ^'^ n—1 ^ 

Then, with these definitions we have the foUowing lemma. 
Lemma 7.2. For any $ G 5 and any t > we have that 

$(Xt,//f) = <I>(Xo,M^)+ f\A^){Xr,fi^)dr+ ^\B<^>){Xr,^i^)dVr+ Ai^'"*) (X, , /i,^)rfr + 

JQ JO Jo 

Moreover, for any T > 0, the following limits hold 

(i~."<i>)(X,,Mf) 



hm I 



dr 

Proof. For p = (p, A) where p = (a, A, cr, f3'~^ ,(3^) G V, define 



= and hm sup E 

N^oo o<t<T 



= 0. 



1 f)"^ f - f)f 

(CVm = 2'^'A^(p) - a(A - A)^(p) 

Then £° is the generator of the idiosyncratic part of the intensity. 
We start by writing that 



k=i •^'■=0 fc=i •^'■=0 

where is a, martingale and 

E I? A^"). ' m").. ■ • • m")..) 



771—1 



AT 

1 

X 



At 



N,3 



1 E {(/:V777)(pf + (/:f^/,.)(p^-)} Mf - 

77=1 

E ^ ((A' A^n. , (/- ^n. ■ • • (/m, mo.) {(^"/- + (^f 

m— 1 

TV 



-^,{{fl,^^^')^,{f„^.^)^...{fM,^^nE)} 



M^. 



At 



Af,4 



^(^*) E £ ((A' /^"). ' (/- A^r), • . . {fM, ),) {(.„(x,)/:f ) J . 

m— 1 



E £ ((/!' A^"). ' (A' Z^").. • • • /^") J {(^f'/™' J • 

m — 1 

To proceed, let's simplify J^ ^. For each / e C°°('P), i > 0, £ N and n e {1, 2, . . . , A^}, define 
(25) JN,n(.t) = -J^ ^(Pt - f [P ■ , \ ' j = „ (£2/, A^t )e - /(Pt )• 

m— 1 
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Then 



< 



7V2 



9A2 



where t ^2l2\ constant from Condition 

Define i(p) =^ A for p = (p, A) G V. Setting 

dip: 



At = Vi(^t) 2^ 

m— 1 
A/ 

= ^i(^*)E 

m— 1 

we have that 



dXr. 



d'~p2 

dXrr, 



1 ^ 



((/i, )^ , (A, A^f )^ . . . (/m, )^) {(Q, A^f )^ (£2/™, A^f )^ - A^f )^} , 



lim E 



r=0 



= lim E 

Af-s-oo 



r=0 



A; 



A; 



= 0. 



Moreover, by Lemma 3.4 in [15] we have that for any T > and any p > 1, there is a constant Co such that 
supo<t<T ^ELi^E < C'o- This and Condition [gj imply that 



lim sup E 

Af->-oo o<t<T 







Collecting things together, we get the statements of the lemma. 
We in particular note the macroscopic effect of the contagion. 

Remark 7.3. The key step in quantifying the coarse-grained effect of contagion was 
average the combination of the jump rate and the exposure to contagion across the pool. 



□ 



Namely, we 



8. Proof of Theorem EH 



Let Qn be the P-law of (X,^^); i.e 



for ah A G ^(i:»Rx£;[0,oo)). Thus Qn G ^(i:»Rx_E[0, 00)) for all iV G N. For w G LIrxbIO, 00), define 

Yt{uj) 1^''a;(t) for all t > 0. 

Also for $ G 5, define the quantity 

rt rt 



(26) 



dcf 



^Yt)-^Yo)- f {A<^){Yr)dr- f {B<^>){Yr)dVr 

Jr=a Jr=0 



Proposition 8.1. We have thatQ^ converges {in the topology 0/ ^(-Drx_e[0, 00))) to the solution Q of the 
(filtered) martingale problem for Af{Y) by i26]) and such that QY^^ — (J^oXttxAo- I'it- particular, Q{Fo = 



cco X TT X Ao} = 1 and for all ^ £ S and < ri < r2 ■ ■ ■ rj ^ s < t < T and C B(R x E), we have 

that Af is a square integrahle martingale with respect to both of the filtrations ^t\jV and ■ Namely, 



(27) 



E 



{Kf{Y)-k1(Y))\{^,{Xr^,^,^^) 



and sup E [Af (F)] < 00 

0<t<T 



Lastly, E[Af (y) V] = 0. 



Proof. The family {/i^jjvGN is relatively compact (as a 15b [0, 00) -valued random variable) by Lemma l7.ll 
Hence {X, /i^jjvGN is also relatively compact. Let {X,fi) be an accumulation point of one of its convergent 
subsequences. Then, {X, fi'^ , J^{B<^){Xr, n^)dVr) will converge in distribution to (AT, fi, {B^){Xr, iir)dVr) . 
This and Lemma [7.2l implv that the process Af (A', /l) will satisfy (|27p. Uniqueness of this martingale problem 
can also be shown as in jl5) . 
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Of course, we also have that for any <J> <E 5, 

EQ[$(yo)] = lim EQ[$(Xo,/x^)] - x tt x Ao) 

which imphcs the claimed initial condition. The rest of the statements are easily seen to be true. 

We next want to identify Q. This will take a couple of steps. 
For notational convenience we shall write 

The next lemma is essential for the characterization of the limit. Its proof is deferred to Appendix 1X1 



□ 



Lemma 8.2. Let W* be a reference Brownian motion and assume that Condition \2.S\ is satisfied. For each 
p = (p, Ao) £ V where p = (a, A, cr, /J*^, there is a unique pair {{Q{t), Xt{p)) : t G [0,T]} taking values in 
K+ X ]R_|- such that 

(28) ^W=/mp.a)6P /3^IEv(Anp)exp[- /* A:(p)dsl|7r(dp)Ao(dA). 



s=0 



s=0 



(29) A*(p) = Ao-a / {XKp) - X)ds + <j ^/X*{^dW: + Q{s)ds + (S'^ X*s{p)dXs. t>0 



s=0 



s=0 



s=0 



(where tt and Ao are as in Condition 



Remark 8.3. For notational convenience we do not write the dependence of Q and A* on X but this is 
always assumed. 

Lemma 8.4. We have that Q ~ 5(x,fi), where for all A G ^{V) and B £ .^(R+), fi is given by 



jltiA X B) [ 



Xi3(At(p))exp 



A:(p)ds 



s=0 



7r(dp)Ao(dA). 



Xa(p)Ev, 

'p=(p,A)e-p 

Proof. Recall ([29]) and the operators £i, £2, -Cf , C% from ([7]) and the definition of Q in ([28 
For any / G C°°{V), 



'p=(p,A)eP 
Using Lemmas IB. II and IB. 21 we obtain 



E- 



/(p. At(p))exp 



A,:(p)ds 



s=0 



7r((ip)Ao(dA). 



p=(p,A)eP 



Ev, 



exp 



A*(p)(is 



7r(rfp)Ao(dA) \ dt 



Ev 



P=(P,A)G-P 



E' 



Vt 



(/:2./)(p,A*(p))g(t)cxp 

(£f*/)(P,A,*(p))exp 



yMds 



s=0 



n{dp)Ao{dX) \ dt 



X*Mds 



s=0 



7r((ip)Ao(fiA) } dVt 



[{C^fjH)E + Q{t) {C2fjH)E + (cf\f,flt)Jdt+(^Cf'f,Pt)^dVt. 



To proceed, define 



dcf 



Git) 

On the one hand, we have that 
Git) = - 



P=(P,A)6-P 



/3^Ev, 



exp 



s=0 



7r((ip)Ao(dA). 



p=(p,A)6-p 



/3^Ev, 



At (p) exp 



Xlip)ds 



s=0 



7r(cfp)Ao(rfA) 
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= (q,A,<7,/3^,/3S) 

On the other hand, by Lemma [8.21 we have 

(30) G{t) = -Q{t). 

Thus, we have that 

d{f,flt)E = { i^if. nt) E + Qii) (^2/, nt)E + (^f /, Mt )^ } rfi + Mt )^ rf^t ■ 

Thus 

Js=0 Js=0 

and hence 5(^x,p.) satisfies the martingale problem generated by A. Of course we also have that p,o ^ it x Ao- 
By uniqueness, the claim follows. □ 



Now we collect our results to prove the law of large numbers given in Theorem 13.11 

Proof of Theorem \3.1\ In Lemma [8.41 we proved that, for any / e C°°(7'), the limiting measure fi satisfies 
the measure evolution equation 



d if, fit)E = {(^1/, f^t)E + (2, (^2/, flt)E + (^f /, f^t)J dt + /, fit)^ dVt. 



From this expression it is immediately derived by integration by parts that, if there exists a solution to the 
nonlinear SPDE then the density of p. should satisfy ([5]). This concludes the proof of the theorem. □ 

Appendix A. Proof of Lemma [HH] 

In this section we prove Lemma 18.21 The proof uses a fixed point theorem argument. For notational 
convenience we sometimes drop the superscript * and p from the notation of (p) and simply write At. 

The square root singularity imposes some technical difficulties in the proof. For this purpose we introduce 
in Subsection lA.ll an auxiliary function t/jjj^x) that will be used later on and study its properties. In Subsection 
IA.2I we study existence and uniqueness and properties of A satisfying (|29p with a given Q{t). 

In Subsection IA.3I we prove the lemma under the additional condition that 60 is bounded. Then, in 
Subsection IA.4I we prove the lemma using Girsanov's theorem for the X process. 

A.l. An auxiliary function. Let < 77 <C 1 and define 

dct 2 Z"''^' ( I 1 dot 

'^''^^^ ^ Inrpi J 1 ' ~Xl-n.v^/^]i^)dz\ dy and g,^{x) = |x|-'0,,(a;) 



-n z 



y=0 IJ z=0 

for all X G M. We note that ipjj is even, so 5,, is also even. Taking derivatives, we have that 

2 1 ■■ 2 1 

for all X > 0. Since g^j = — < 0, 5^ is non-increasing. For x > y/rj, 

^ In 77^/^ - In?? 
g,j{x) = 1 - 2 j-^ = 0, 

so in fact is nonncgative on (0, cxj) and it vanishes on [y^, 00). Thus 5^ is nondecreasing and reaches its 
maximum at y^. Since g,j{0) = 0, we in fact have that 

< < .9,, (07) 

for all .T > 0. Since g,j is non- increasing on (0, 00) and gri{x) = 1 for .t e (0, 77), we have that grj{x) < 1 for all 
X e (0, ^/rj), so grii\/v) ^ Since g,, is even, we in fact must have that |g,,(x)| < y/rj for all x G M. Hence 

\x\ < i^rjix) + y/rj 
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for all a; G R. We finally note that 



(31) 



< 



In 77 



1 1 

|a;| ' 77 



and that xijjjjix) > for all x G 



A. 2. The uncoupled linear case. Let ^ be a {^^f}t>o-predictable, nondccreasing, bounded and right- 
continuous process such that ^0 = 0. Consider the SDE 



(32) 



dXt = -a{Xt - X)dt + a^/Xt^OdWi + d^t + P'^XtdXt t > 
Ao = Ao- 

Then, as in Lemmas 3.1 and 3.2 of |15] . we get that (|32p has a unique nonnegative solution such that 
suPtg[o.T] lEilAtl''] < 00 for all r > and q>l. 

A. 3. Proof assuming that boi^x) is bounded. In this subsection we prove Lemma \S72\ assuming that bo{x) 
is bounded. The proof uses a fixed point theorem argument. 

The main condition of this subsection is that there exists a A/ < 00 such that 



(33) 



sup 1 60 (x) I < M 



Let q> 1 and S''(]R_|-) be the set of valued, adapted, continuous processes {Xt}t£[o.T] such that 



|A|| 



sup E|At|« 
a<t<T 



1/9 



< 00 



The space S''^(]R-|_) endowed with the norm is a Banach space. 

Let us consider now a nonnegative process Ut{p) € S'''(R+). Set 



(34) 



p=(p,A)e-p 

p=(a,A,a,/3^,/3S) 



V i cxp 



Us{p)ds 



7r(dp)Ao(dA). 



and consider ([5^ with ^ in place of ^{U). We are going to prove that the map $ defined by A = ^{U) 
through with U £ S^(R+) is a contraction on S'^(]R+) equipped with the norm 

l|A||f,i < 00 

locally in t which can also be extended to any arbitrary T. 

We collect in the following lemma some important properties of ^ as defined though (|34p . 



Lemma A.l. Let a nonnegative process Ut{p) G S'*(M+) be given and let K be the constant from Condition 
2.^ The map ^{U)t, as a function oft, is continuous, non- decreasing, positive, bounded uniformly m i G M_|_ 
by K and satisfies £,{U)o = 0. 



Proof. All the statements are obvious. 



□ 



Fix [/((p), U[{p) G S'^(]R+) and write for notational convenience = ^(C^)t and = £,{U')t- The process 



Zt — Xt — X[ satisfies 



Zt 



ZAs 



+ T / ( VaT - v/aT) dWs + / (Cs - dQ + /? 



ZAX, 



Next, we apply Ito formula to iprj{x) defined in Subsection I A. II with x = Zt, getting 



i^riiZt) 



ZMZs)ds+ / i;:JZs){d^s-dQ 



(35) 



+ Y (v^-VaT) ^y^iZs)ds + a (VaT- v/Ar)v;(^s)dW^s 

+/3^ / Po{Xs)Zs^'^{Zs)ds + p'' [ <T,{X,)Zsij'^{Zs)dV,+ [ {p^(T,(Xs)Zsfi^':^{Zs)ds 



For notational convenience we define the quantity AjK = Yj — for any given couple of stochastic 
processes Y, Y' . 

We have the following three technical Lemmas IA.2IIA.4l 

Lemma A. 2. For any q > I and t > we have 

E|AtC|« < A-J /" E /" |A,t/(p)|«ds7r(dp)Ao(dA) <tA-? /" ||A.;7(p)||«^^(dp)Ao(dA), 

"'p=(p,A) Jo Jp={p,X) 

where K is the constant from Condition 

Proof. By the definition of £,{U) from (|34p we get 



P=(P,A) 



13^ [Ev< exp 



U'Mds 



s=0 



exp 



Us{p)ds 



s=0 



7r(dp)Ao((iA). 



Recall the trivial inequality e" — e < maxje"^, e }\a — 6|. Since C/s(p), [/^(p) > 0, we get that 



exp 



u'Mds 



=0 



exp 



Us{p)ds 



s=0 



< 



\Us{p)-U:{p)\ds 



The latter, Holder inequality and Fubini's Theorem imply 

ElAt^l'^ < K" f [ E|A,[/(p)|%(dp)Ao(dA)d 

< tRi 



sup E|A,[/(p)|%(dp)Ao(dA) 

p=(p,A) 0<s<t 



This concludes the proof of the lemma. 
Lemma A. 3. The following bound holds 

(36) EiJjrjiZt) < Ci f E\Zs\ds + C2 



□ 



||A.C/(p)||,,^(dp)Ao(dA) + C3 

'P=(P,A) 

where Ci = Ci {K, M), C2 = C2{K, t) and C3{K, t, r]) ^0 as rj ^ 0. Here K is the constant from Condition 
\2.<A and M as in 



Proof. The proof of this lemma follows by bounding each term on the right hand side of psp separately using 
the bound from Lemma [A . 2 1 and the bounds for the first and second order derivatives of iprj, i-e. < 1 

and ()3ip respectively. In particular we have the following. 
Taking expected value in psp we obtain 







aE Zs^'^{Zs)ds + E pHo{Xs)Z,^[^{Zs)ds + E ^'^{Z s) {dt, ^ dQ 



(37) 



(TaT-V^) i^l'j{Zs)ds + E {l3'^ao{Xs)Zs)^4',iZs)d. 



Let us now bound each term on the right hand side of ([37| . 

Due to the boundedness condition on bg we have for the first and second term 



aE Z,i''^{Zs)ds + E p^bo{Xs)Zs4>'^{Z,)ds < {a + \P^\M) E\Z,\ds 



(38) 

Regarding the third term we recall the properties of ^ outlined in Lemma lA.ll and that |7/'^(x)| < 1 for 
all 77 > and x G M. Then, approximating by simple processes, we get that that there exists a constant Cq 
such that 

<Co sup EiCs-e^i 

0<s<t 



^/(Z.)(d6-0 
The latter display and Lemma lA.21 implv that 



(39) 



< CotK 
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P=(P,A) 



|A.[/(p)||,i^(dp)Ao(dA) 



Relation ([31]) gives for the fourth term 



(40) 



-E 



< 



(E|Z,|V;'(Z.)) ds< 



In rj 



Relation ([3T|) and the condition E J* {ao{Xs))'^ ds < oo give for the fifth term 



(41) 



E 



(/3^ao(X,)Z,)'v4(^s)rf^ 



< 



In 77^^ 



Collecting now terms (|55| - (|4ip we obtain the statement of the lemma with 
Ci = a+|/3^|M 
C2 = CotK 



□ 



Lemma A. 4. For any t > 0, there exists a positive C{t) such that 



||A.A||, i <C(t) / ||A.C/(p)||,,i7r(dp)Ao(dA). 

Moreover, C(t) is continuous, increasing in t and satisfies limt^oC'(0 ~ 0- 

Proof. Applying |a;| < ipnix) + y/rj to x = Zt and using Gronwall's Lemma we obtain 



< 



Co 



P=(P,A) 
t 







||A.(7(p)||,_i7r(dp)Ao(dA)+C3(t,77)+7yi/2+Ci e<^^(*-^) (Cs{s,v) + v'^^) 
C2\\A.Uip)\\s.i<dp)KidX)ds 



ds 



P=(P,A) 



where the constants Ci , C2 , C3 are as in Lemma IA.3I Taking r/ J, we obtain 

nZt\<C{t) f ||A.(7(p)||,i7r(dp)Ao(dA) 

Jp=(p,A) 

where C{t) — Co[l + e^^* — e^^^] t — Coe^^*^t for some constant Cq > 0. This concludes the proof of the 
lemma. □ 



We collect the previous results for the proof of Lemma 18.21 in the case of a bounded drift . 
Lemma A. 5. IfbQ{x) satisfies LS3\) . then the statement of Lemma \8.2\ is true. 

Proof. Using Lemma lA.4| a standard Picard iteration procedure shows that there exists a fixed point A* of 
$, i.e. \l = $t(A*), for t e [0,ti] such that C{ti) < 1. 

Let us show that this fixed point is necessarily unique. Indeed, suppose that A*, A are both fixed points. 
Then, Lemma lA.41 implies that 



A*(p)-A(p) <C[t) A*(p)-A(p) ^(dp)Ao(dA) 

*'l ■/p=(p,A) 

By integrating, and using the condition C{t) <\ for t e [0,ii], we immediately obtain that 

A*(p)-A'(p) ^(dp)Ao(dA) =0. 



P=(P,A) 

Therefore, for every t € [0,ii] and p e V we should have A^ (p) = Aj(p) almost surely. This gives uniqueness. 

For the general case, we only have to subdivide the interval [0,r] into a finite number of small intervals. 
Thus, there is a unique fixed point to the equation A^ = $t(A*) for all t e [0,T]. Notice that the equation 
Aj = $t(A*) can be trivially written as the pair of the coupled equations (|28|) - (l29l) . Clearly (5(t) is nonnegative 



since A* is nonnegative. This concludes the proof of the lemma. 
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□ 



A. 4. Proof of Lemma l8.2l In this section we prove Lemma lS^ By Lemma lA.SI we know that the statement 
of Lemnia l8.2l is true if bo{x) is bounded. Our strategy is to apply this resuh for the case 60(2;) = and then 
to generalize using Girsanov's Theorem. 

Recall the map $ defined by At — $t(A) through p2p - ([M|) . For notational convenience we shall write 
A° for A( = 4>t(A°) to emphasize the fact that we are considering with bo{x) = 0. We keep the 

notation A for the case of a general bo{x). 

Let u{x) be such that <7o{x)u{x) = ^bo{x) and define the quantity 

For given Aq G S'^(M+) we define iteratively the sequence 

A„+i,t = $t(A„) 

We shall prove that the sequence {A„}„gN is a Cauchy sequence in probability uniformly in t € [0,T]. In 
particular we have the following lemma. 

Lemma A. 6. Assume that there is a p > 1 such that EM^ < 00 and let q > 1 be such that l/p + 1/q = 1. 
For every e > 0, there exists N < 00 such that for all N < n < m 

sup ¥[\XT-X'l\>6]<Co{e/6f' 

0<t<T 

for every 6 > and for some constant Cq > independent ofn,m. 

Proof. As we saw in the proof of Lemma [A. 51 it is enough to consider the equation At ~ $t(A). Let A*, A 
be two solutions of the equation At ~ ^t{X)- 

By Girsanov's Theorem on the absolute continuous change of measure in the space of trajectories, Holder 
and Chebychcv inequality, we have the following 



sup P[|Ar-Ar| ><5] - sup E[x{|Ar-A?|>5}] 

0<t<T 0<t<T 



sup E 

0<t<T 



X{|A-"-A-"|>5} 



Mt 



< sup P [|A°''" - A°'"| > 6] {EMP f^P 

\0<t<T J 

< |^^^po<.<TiE[K-"'-Ar"|]y^^^^^^^^,,, 

where A°'™,A°'" satisfy A°'™ = «>t(A°'™) and Aj'" = $t(A°'"). 

Then Lemma lA.41 implies that supQ<t<TE [|At'™ — At'"|] can be made arbitrarily small by choosing n,m 
large enough. This together with the boundedness of EM^ conclude the proof of the lemma. □ 



Now we are in position to prove Lemma 18.21 

Proof of Lemma \8.2l By Lemma IA.6I we know that there is a nonnegative solution to the equation At = 
<f>t(A). Uniqueness is shown in a similar fashion. This concludes the proof of the Lemma. □ 

Appendix B. Some auxiliary results 

In this section we prove some auxiliary lemmata that were needed in various places in the paper. These 
results arc interesting themselves in their own merit. 

Lemma B.l. IfY is an integrable ^t — n^easurahle random variable, then E \Y\Vt\ = E \Y\V]. 
Proof. Notice that V = Vt V where 

V+ ^a{Vt+s-Vus>0). 

Since V is an ^t— adapted Brownian motion, we get that C V is independent of J^f Therefore, we have 

E[y|Vt] = E y|vtYVt+] =E[y|v]. 



□ 



Lemma B.2. Consider a cadlag and adapted process {Yut> 0} such that E /p Y^dt\ < oo for every 
T > 0. Then we have 



E 



Y.^dV, 



Vt 



= / E 



V,, 



dK. 



Proof. Given an arbitrary u S L°°([0,i];R), i.e. supg<j,<( \us\ < oo consider the quantity 
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It is easy to see that Zt satisfies the SDE 



Zt = l 



Therefore, on the one hand we have 

ft 



ZtE 



YdV, 





















Vt 




= E 


Zt f Y.dVs 




= E 


V 








Jo 






. j( 
















= E 


i / ZsUsE 


Ys 




ds 






. Jo ^ 









(42) 

On the other hand 



E 













r /•* 








Zt f E 


Y, 




dVs 


= E 


i / ZsUsE 


Ys 




ds 


Jo 



















Thus we arrive at the equahty 



Vt 



Zt 



Ys 



Vs 



dVs 



YsdVs 

Since, this statement is true for every u e L°°([0,t];M) we conclude the proof of the lemma. 
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